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Abstract — The influence of interface roughness scattering (IRS) 
on the performances of silicon nanowire field-effect transistors 
(NWFETs) is numerically investigated using a full 3D quantum 
transport simulator based on the atomistic sp 3 d s* tight-binding 
model. The interface between the silicon and the silicon dioxide 
layers is generated in a real-space atomistic representation using 
an experimentally derived autocovariance function (ACVF). The 
oxide layer is modeled in the virtual crystal approximation (VCA) 
using fictitious SiC>2 atoms. (110) -oriented nanowires with dif- 
ferent diameters and randomly generated surface configurations 
are studied. The experimentally observed ON-current and the 
threshold voltage is quantitatively captured by the simulation 
model. The mobility reduction due to IRS is studied through 
a qualitative comparison of the simulation results with the 
experimental results. 

Index Terms — Si gate-all-around nanowire transistors, atom- 
istic, full-band simulations, interface roughness scattering 



I. Introduction 

AS the dimensions of the planar metal-oxide- 
semiconductor field effect transistor (MOSFET) 
are approaching the nanometer scale, their performances 
become limited by the short channel effects (SCE), an 
increasing OFF-state current, and a poor electrostatic control 
of the channel through a single-gate contact. 

The nanowire (NW) field-effect transistor has been identi- 
fied as a promising candidate to overcome these issues and 
build the next generation switch [1], [2]. Several experimental 
studies have shown that NWFETs have excellent electrical 
characteristics [3]-[5], especially in a gate-all-around (GAA) 
configuration that provides a superior channel electrostatic 
control [6], [7]. Among other advantages, GAA NWFETs 
exhibit the smallest SCE as compared to other multi-gate 
(MG) structures [8], [9]. Despite the better electron transport 
properties of III-V materials, Si has remained the material of 
choice to fabricate GAA NWFETs due to its well-understood 
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characteristics and its compatibility to the conventional CMOS 
technology [10], [11]. 

In line with experimental efforts, simulation and model- 
ing approaches have been developed to help design better 
performing NWFETs and understand their physical behavior. 
Until very recently, drift-diffusion (DD) and Monte Carlo 
(MC) approaches have been successfully used to simulate large 
devices. However, at the nanometer scale it is necessary to go 
beyond DD and MC models and use a quantum mechanical 
approach that captures the wave nature of electrons [12]. 

Nano-scale devices are characterized by 3D material varia- 
tions at the atomistic scale so that an atomistic description 
of the simulation domain has become indispensable. Non- 
parabolicities and anisotropics of the bandstructure strongly 
affects the transport characteristics of electrons in NW tran- 
sistors [13], [14]. Furthermore, the strong confinement of elec- 
trons in nano-devices lifts degeneracies due to valley splitting 
[15]. These effects cannot be captured by the effective mass 
approximation. Therefore, it is inevitable to go beyond the 
effective mass approximation and to use a full-band approach 
such as the nearest-neighbor empirical tight-binding (ETB) 
model [16]. 

With the recent development of OMEN [17], a 3D full -band 
quantum transport simulator based on the atomistic sp 3 d 5 s* 
TB model [18], it has been possible to explore nano-scale 
NWFETs. However, the computational burden of such an 
approach is so intense that large computing resources and an 
efficient parallelization scheme of the work load is required 
[19]. 

The focus of this paper is on interface roughness scattering 
(IRS) that is, in conjunction with electron-phonon [20] and 
impurity scattering [21], crucial to understand the electrical 
characteristics of ultra-scaled NWFETs [22]. At the nanometer 
scale, the details of the NW surface become more important 
[23], [24]. In the ballistic transport limit, the normalized ON- 
current of NWFETs is expected to increase as the diameter 
decreases. Experimental data shows that the ON-current in- 
creases with decreasing diameter until 3 nm [22]. However, 
decreasing the diameter below 3 nm reverses the trend and re- 
duces the current density. This behavior has been attributed to 
electron-phonon scattering and interface roughness scattering 
[22]. Here we examine this experimental reasoning through 
modeling and simulations. 

In traditional MOSFETs, IRS has been identified as one 
of the most important scattering mechanisms [25]. Because of 
the strong confinement of carriers in inversion layers that push 
them close to the semiconductor-oxide interface, the effective 
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mobility of MOSFETs is limited by IRS at high electric fields 
or at high electron densities [26], [27]. 

IRS in classical MOS type devices is treated via perturba- 
tion theory [28]. It is represented by an IR-limited mobility 
depending on the root-mean-square (RMS) of the roughness 
amplitude and its correlation length [29]. Early 1-D quantum 
device simulations of IRS in the NEGF formalism also used 
a perturbation treatment [30]. In 3D nanowire devices, an 
explicit IR representation is feasible and could represent the 
device- to- device fluctuation. Recent theoretical studies have 
shown that mode-mixing due to IRS is reduced in NWFETs 
compared to conventional MOSFETs [31 ]— [33] and the po- 
tential fluctuation throughout the NW dominates IRS in small 
cross-section NWs [33]. However, these studies have been 
limited to a continuum representation of rough surfaces with 
(100) as transport direction and cross-sections larger than 3 x 
3 nm 2 . 

Atomistic full-band simulations offer a more realistic repre- 
sentation of rough surfaces for any crystal directions. Atom- 
istic simulations show that the ON-currents and the mobilities 
of NWFETs are reduced significantly by IRS when the di- 
ameter of NWs is below 3 nm [34], [35]. These simulations 
have been so far limited to ultra-narrow NWs without any 
representation of the SiC>2 layer, i.e. hard wall boundary 
conditions are applied to the silicon surface and the electrons 
do not penetrate into the oxide layer. 

In this work, IRS is investigated in GAA <110> NWFETs 
with different diameters (2 nm - 4 nm) through an atomistic 
representation of the rough surface. Contrary to previous 
studies, the S1O2 layers are modeled explicitly and taken into 
account not only in the Poisson equation, as in Ref. [36], but 
also in the quantum transport calculation. The oxide layers 
are modeled in the virtual crystal approximation (VCA) where 
fictitious SiC>2 atoms are used. Disorder is introduced through 
a disordered spatial representation of Si and "SiCV atoms. 
The inclusion of the wavefunction penetration into Si02 allows 
us to investigate the ON-current and mobility reduction due 
to IRS in NWFETs and to compare simulation results to the 
experimental work in Ref. [22]. 

The paper is organized as follows: In Section II, the 
simulation approach is introduced with an emphasis on the 
generation of rough surfaces in circular nanowires. In Section 
III, the simulation results are presented and discussed. 

II. Method 

A. Interface roughness model 

The statistical properties of the interface between Si and 
SiC>2 is characterized by the autocovariance function (ACVF) 

where A m is the root mean square (RMS) of the interface 
height and L m the correlation length of the rough surface. 

Here, s(r^) is a random number representing the amplitude 
of the rough surface at position T\j where a Si atom is located 
and connected to SiC<2 atoms as shown in Fig. 1(a). The 
quantity s(ri) can be defined as the distance between the 
position of the Si atom and the reference plane. This model 
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Fig. 1 . (a) Atomistic representation of the cross-section of a perfect nanowire 
along the (110) direction, (ref.: reference plane, Cy/Cy : ACVF for Si/SiC>2, 
rilr 1 ^. position vector of i-th Si/SiC>2 atom) (b) 2D surface profile generated 
for a nanowire of diameter 3 nm with Z\m — 0.14 nm, Lm — 0.7 nm 
for the NWFET structure presented in Fig. 3 (y' is a coordinate along the 
circumference direction.) 

has been developed from an experimental study [37] where 
the measurement on the roughness amplitude is performed at 
an atomistic level. It is used in most of the simulation studies 
on IRS in NWFETs [31]-[33], [36]. Notice that the y/2 term 
in the numerator of the exponent in Eq. (1), which is not 
present in Ref. [36], is included to correctly fit the model to 
the experiment as in Ref. [37]. 

From Eq. (1), one can calculate s(rj) for each Si atom at 
the interface following the procedure in Ref. [38]. Then, s(r^) 
for the Si02 atoms at the interface should be calculated by 
averaging the s(ri) of Si atoms that are connected to a Si02 
atom located at r\ as shown in Fig. 1(a). This introduces an 
error into the ACVF. For example, for a correlation length 
L m = 0.7 nm, the error due to this approximation is ap- 
proximately 1.88% only (see Appendix A for a more detailed 
derivation). 

After s(r ?; ) and s(r£) are calculated for the atoms located 
at the first Si-SiC<2 interface layer as shown in Fig. 1(a), the 
s(-) for atoms at the second layer are calculated in the same 
manner. Then, Si atoms are replaced by SiC<2 atoms and vice 
versa depending on the value of s(-) and a selection criterion. 
The criterion is described as follows: let d(ri) be the distance 
from the atom considered to the reference plane located in 
between Si and S1O2 atoms. Then, if s(ri) > d(r;), the Si 
atom at should be replaced by S1O2 and if s(r^) < — d(rQ, 
the Si02 atom at r\ should be replaced by a Si atom. 

The rough interface of a NW of diameter 3 nm and length 
15 nm is projected to a two dimensional (2D) plane and plotted 
in Fig. 1(b). The rapid short-range fluctuation of the surface 
profile is the characteristic of the chosen exponential model 
rather than a Gaussian model [37]. 

At this point, it is necessary to verify that the rough interface 
profile generated with the algorithm explained above has the 
same statistical properties as the theoretical ACVF (TACVF) 
presented in Eq. (1). It is not clear in the literature [31]- 
[33], [39] whether the generated rough interfaces are correct 
representations of theoretical ACVF or not. Only in Ref. 
[40], the one dimensional (ID) ACVF of a single sample is 
compared to its theoretical value and it is shown that the first 
few coefficients of the ACVF from the lowest order agree with 
the TACVF. 

In this paper, 2D ACVFs of the generated rough interface 
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Fig. 2. (a) The average of ACVFs for 100 samples to x direction, and (b) 
to the circumference direction (y') are displayed. The halves of the graphs 
(a,b) to the negative direction are omitted because of their symmetry. All the 
calculations are done for NWs with diameter 3 nm and gate length 15 nm 
(see Fig. 3). 

are compared to the TACVF after averaging ACVFs of 100 
samples. The ACVFs can be obtained either from s(r^) 
(labeled ACVF1) or from the position of the rough interface 
atoms (labeled ACVF2). ACVF1 should agree with TACVF 
and ACVF2 with ACVF1 and consequently with TACVF. The 
following equation is used for 2D ACVF calculation [41] 

N m -i-iN y -j-i 

Qj = AT AT y ' y (fm.n ~ f)(fm+i,n+j ~ /)> 
y m— n— 

(2) 

where f m , n is the amplitude of the rough surface perturbation 
projected onto the regular grids, and / is the average value of 
fm.n- N x and N y are the numbers of grid points on the x and 
y' axis (see Fig. 3 for a reference), respectively. 

Since it is not easy to directly compare 2D ACVFs with each 
other, cutting sections of 2D ACVFs through the origin to x 
and y' direction are compared in Fig. 2. The ACVF1 and the 
ACVF2 agree with TACVF for both x and y' direction when 
they are averaged over the entire rough NW samples. This 
validates the algorithm used for generation of rough interfaces 
in this work. The importance of the correct representation of 
interface roughness can be illustrated in the effect of the RMS 
value (=^/ACVF(0, 0)) on the drain current from Fig. 11(d). 

B. S1O2 model 

Previous attempts to model oxide materials in the tight- 
binding framework have been made by assuming a well- 
defined crystallographic structure of SiC<2 such as /3-quartz, 
tridymite [42], or /?-cristobalite [43]. All these approaches 
have in common a sp 3 description of SiC<2 up to second 
nearest-neighbor interactions for the oxygen atoms. Their 
application is limited to one-dimensional Si-SiC<2-Si structures 
with a semiconductor-oxide interface perpendicular to the 
<100> crystal axis only. 

As an alternative and to go beyond one-dimensional <100> 
transport we introduce a simplified, but multi-functional model 
to describe silicon dioxide. Fictitious Si02 atoms are consid- 
ered and treated in the virtual crystal approximation. They are 




Fig. 3. The simulated structure of NWFET: A circular NW with varying 
diameter (Dsi=2, 2.5, 3, and 4 nm) is used. The default values of the other 
parameters are displayed in the figure. 

arranged as a diamond crystal with the same lattice constant 
ao = 0.543 nm as the Si atoms. A nearest-neighbor sp 3 tight- 
binding model is used to describe the electrical properties of 
the Si02- The band gap and conduction band offset are fitted 
to the experimental values with the relations E g — 8.9 eV and 
AEcb = 3.15 eV with respect to Si [44] as well as an electron 
effective mass m e = 0.44mo [45]. A detailed description of 
the model is presented in Appendix B. 

III. Results 

The NWFET structure considered in this work is illustrated 
in Fig. 3. A m = 0.14 nm and L m = 0.7 nm are used as default 
parameters which are experimentally measured in Ref. [37] for 
a 2D Si/Si02 interface. Crystal orientation (110) is selected 
as the transport direction because of its superior transport 
characteristics compared to other crystal orientations [13] and 
the highest immunity to interface roughness scattering [36]. 
Two important parameters the diameter of the NW channel 
Dsi (depicted in Fig. 3) and A m are varied to study their 
impact on the NW performances. The diameter of silicon NW 
channel is varied from 2, 2.5, 3, to 4 nm and the RMS of 
the roughness amplitude is varied from 0.14, 0.2, to 0.3 nm. 
Strain is not considered in this work. The calculated Id-Vg 
characteristics of perfect and rough NWFETs of 100 samples 
with a diameter of 2 nm are shown in Fig. 4. 

Fig. 4 presents two effects of IRS on Id-Vg- The first effect 
is the positive shift of the threshold voltage and the second 
effect is the reduction of the ON-current. The positive shift 
of threshold voltage is ascribed to the fact that the diameter 
of rough NWs throughout the channel is reduced (on average) 
from the perfect soft wall NW due to rough interfaces [46] 
and that the density of states (DOS) in the energy spectrum 
is raised as discussed in Ref. [31]. The amounts of threshold 
voltage shift on average are 49.2, 12.5, 8.1, and 10.0 mV for 
rough NWs with diameter 2, 2.5, 3, and 4 nm. The standard 
deviation a of the threshold voltage shift is 13.6, 10.1, 5.5, 
and 2.0 mV for each diameter from 2 to 4 nm. 

The impact of IRS on the ON-current is measured by the 
reduction of the ON-current (/ pe rf — ^scatt) /-fperf x 100 where 
iperf is the drain current for perfect NWs with softwall BC and 
^scatt is the drain current for rough NWs. The ON-currents of 
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Fig. 4. Id vs Vq for NWFETs of Dsi = 2 nm and total diameter Dj = 
X?Si + *Si0 2 x 2 = 4 nm. I cr it(= Dsi x 10 — 7 A) is the critical current for the 
threshold voltage calculation. The gate work function for this simulation is 
chosen to be 3.95 eV. Vb is set to 0.6 V. The error bars indicate the standard 
deviation of the current in each bias point. 



rough NWs are reduced from perfect NWs by 20.5 % (Z?si = 
2 nm), 16.5 % (D si = 2.5 nm), 9.0 % (D si = 3 nm), and 3.8 
% (Dsi = 4 nm) due to IRS on average. 

The effect of hardwall/softwall BC on the threshold voltage 
is noticeable as shown in Fig. 4. It can be explained from 
Fig. 5 where the electron density is plotted through the center 
of the hardwall/softwall perfect NWs at (z — 0, x — Lq/2). 
An increase of the carrier concentration near the interface 
between Si and SiC>2 layer in the softwall NW reduces 
the effective oxide thickness leading to the decrease of the 
threshold voltage. 

The standard deviation of the threshold voltage due to IRS 
in the hardwall rough NWs (for 100 samples) for diameter 2 
nm is calculated to be 14.3 mV, that is, 5% increase compared 
to that in the softwall rough NWs. The effect of BC on the 
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Fig. 5. Electron density along j/-axis in the middle of the channel at 
(z = 0, x = Lq/2) of a nanowire transistor with hardwall or softwall BC in 
the Si and Si02 interfaces. The inset describes the Si-SiC>2 interfaces around 
y = 1 nm. 
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Fig. 6. (a) The threshold voltage and (b) the ON-current at Vg — Vth = 0.4V 
vs diameter for the perfect and rough NWFETs. The experimental ON-current 
and the source series resistance Rs are extracted from Ref. [22]. The results 
with phonon scattering are extracted from Ref. [20]. The inset in (a) magnifies 
the threshold voltage trend in the square box. 



ON-current is even larger. The ON-current reduction in the 
hardwall NWs is 30% which is 9.5% higher compared to 
20.5% reduction in the softwall NWs. The standard deviation 
of the ON-current in the hardwall NWs is 0.179 mA/pm 
compared to 0.141 mA//im in the softwall NWs, that is, 26.7% 
increase. As the effective cross-section of the hardwall NWs 
is smaller than that of the softwall NWs, the fluctuation of the 
conduction band edge due to IRS is larger in the hardwall NWs 
and hence the calculation using the hardwall BC overestimates 
the effects of IRS. 

The threshold voltage and the ON-current for NWFETs with 
different diameters are plotted in Fig. 6. The threshold voltage 
is calculated using the constant current method [31], [47] with 
a critical current I C rit(= D si x 10~ 7 A). The threshold voltage 
shift due to phonon scattering is extracted from Ref. [20]. Be- 
cause phonon scattering reduces the drain current for the whole 
range of the gate bias, the threshold voltage is also shifted 
to achieve the same critical current. The threshold voltages 
calculated for the perfect NW and the rough NWs have similar 
values and slopes when the diameter of the NWs is larger than 
2.5 nm, and show similar trends as the experimental values 
extracted from Ref. [22]. However, as the diameter decreases 
below 2.5 nm, the slope of the Vih--Dsi curve for the perfect 
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Fig. 7. (Left) the current spectra (transmission X (/l — /r)) with the Fermi 
level in the left contact _E/l and (right) the density of states in a log scale 
across a rough nanowire FET resolved in energy E and longitudinal coordinate 
x at the ON-state (Vq — Vth ~ 0.4 V) overlapped with the conduction band 
edge E c - 



NW deviates from the experimental values. This implies that 
to get a consistent trend or quantitatively matched threshold 
voltage, one needs to include IRS in the transport calculation 
when the diameter of NW is scaled below 2.5 nm. 

The drain current calculated at Vg — V t h = 0.4V, i.e., the 
ON-current is shown in Fig. 6(b). The ON-current reduction 
due to IRS at a diameter smaller than 3 nm is significant 
and partly explains the decrease of the ON-current in the 
experimental results [22]. 

The ON-current with phonon scattering is calculated from 
the ballisticity in Ref. [20] of a very similar structure where 
only the doping density in the source/drain was different. 
In Ref. [20], the ballistic current is calculated using the 
electrostatic potential extracted from the phonon scattering 
simulation to exclude the contribution of the phonon scattering 
in the source/drain extension region. 

Since fully atomistic NEGF calculations are extremely com- 
putationally intensive, requiring thousands of CPUs for a sin- 
gle I-V, we cannot afford to perform statistical sampling with 
IRS. We therefore extract the effect of phonon scattering via 
the ballisticity obtained in Ref. [20]. However, the simulations 
in Ref. [20] are done for (100) -oriented NWs whereby our 
simulations are conducted for (llO)-oriented NWs. From the 
fact that (llO)-NWs have approximately 20% larger ballistic- 
ity than (100)-NWs when phonon scattering is included [48], 
the ballisticities of (llO)-NWs are approximated from those of 
(100)-NWs. The relationship 1/T IR+Phn = l/T IR +l/T Phn -l is 
used where the transmission T is calculated from the ballistic- 
ity factor 13 using the relation T mvi m = 2/?m/phn/(l + /Win)- 

The ON-current reduction due to phonon scattering exceeds 
that of IRS, and the drain current including both phonon and 
IR scattering (with series resistance) are in accordance with the 
experimental results. Therefore, phonon scattering should also 
be included in the transport calculation to get a quantitative 
prediction for the ON-state characteristics of NWFETs. 

Fig. 7 conveys detailed information about the effect of IRS 
on the current spectra and the DOS. IRS causes reflection of 
electrons and reduces the current spectra in the energy range 
where the peak is located in the perfect NW. A localization 
effect of the DOS, indicated by two circles, due to the 
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Fig. 8. (Top) the electron density profile through the rough NW with 
higest/lowest drain current, (bottom) transmission through the perfect NW, 
the rough NW with the highest/lowest drain current at Vq = Vb = 0.6V, 
and (inset) Id — Vq plot. 



fluctuation of the conduction band edge Ec in the channel 
causes this reduction of the current [34], [46]. The current 
spectra and the DOS plot also shows that the relevant energy 
region of the electron transport includes subbands that are 
higher than the second subband. This supports that a full- 
band model is necessary to correctly understand and model 
these devices. 

The transmission curves for the perfect NW, the rough NW 
with the lowest/highest drain current in Fig. 8 shed light on 
the impact of IRS on the electron transport characteristics of 
the nanowire transistor. The decrease in transmission due to 
IRS throughout the whole range of energy is observed. For 
the rough NW with the smallest drain current shows more 
localization effect in real space compared to the rough NW 
with the highest drain current, as can be seen from the top 
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Fig. 9. The ballistic efficiency /3 sal which is the ballistic ON-current divided 
by the average ON-current with IRS, the phonon scattering (extracted from 
Ref. [20]), or both the IR and the phonon scattering included. 
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Fig. 10. (a) The IR-limited mobility ^ir is calculated for the electron density 
N a = 3.5 X 10 6 cm _1 from the mean of the effective mobility including IRS 
Mcff.lR an d the ballistic mobility through the Matthiessen's rule fi^ = 
^efflR — ^ball ■ C 5 ) ^ e ballistic and effective mobility calculated for gate length 
Lq = 30nm are compared to the experimental results which are extracted 
from Ref. [22]. Effective channel length is taken into account by reducing 
the channel length below the gate length according to the experimental result 
[22]. 



figure in Fig. 8. 

To further understand the importance of IR, the ballisticity 
is calculated for rough NWs with D^i = 2nm after A m is 
increased to 0.2 nm and 0.3 nm. It is then compared to the 
ballisticity including phonon scattering extracted from Ref. 
[20]. As shown in Fig. 9, the impact of IR on the ballisticity is 
smaller than that of phonon scattering for all diameters, but it 
becomes comparable to the phonon effects as A m increases 
to 0.2 nm and starts to exceed the phonon effects when 
A m = 0.3 nm. 

As the ON-state characteristics of a transistor are also 
related to their effective mobility [49], [50], the electron 
mobility reduction due to IRS at low drain bias should be 
understood in relation with the ON-current reduction at high 
drain bias. The low-field effective mobility in NWFETs is 
calculated from the electron density and the conductance using 
the expression [32], [33] 

Glin^ch 

Mcfi = — T- T — (3) 

where Gy m is the conductance in the linear region i.e. at low 
drain bias (Vu = 5 mV in this work) and iV c h is calculated 
by integrating the electron density in the subsection of the 
channel under the gate where the electron density is nearly 
uniform [33] as illustrated in Fig. 11(c). 

Assuming that the ballistic mobility is the effective mobility 
for a perfect NW, the IR-limited mobility /iiR can be calculated 
from the Matthiessen's rule, namely, fj,^ = Mini - Mbaii wnere 
MeffjR is the effective mobility for rough NWs. The IR-limited 
mobility /iir plotted in Fig. 10(a) is calculated from the mean 
of the effective mobility. It decreases as the diameter of the 
NW decreases. This is expected because IRS increases as the 
diameter decreases since the electrostatic potential fluctuates 
more in NWs with smaller diameter [51] and it is consistent 
with the ON-current behavior in Fig. 6(b). The rate at which 
the average of ^ decreases is approximately which is 
much smaller than for small electron density as described 
in Ref. [51]. 

To calculate the phonon-limited mobility, the resistance of 
the nanowire transistor is calculated for the same structure 
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0.3 



Fig. 1 1 . (a) The IR-limited mobility /jjr compared with the phonon-limited 
mobility ^ p hn> (b) the total effective mobility including only IR or both the IR 
and the phonon scattering through Matthiessen's rule compared to the ballistic 
mobility rtiall. ( c ) ID electron density along the channel with the vertical bars 
defining the subsection of the channel for mobility calculation, and (d) the IR- 
limited ON-current compared with the ballistic current are displayed. All the 
figures are plotted for the nanowires with diameter D$i=2 nm. The electron 
density is set to 7V S = 3.5 X 10 6 cm _1 except for (c) where the gate/drain 
voltage is fixed. 



as shown in Fig. 3, but for several different channel length 
of nanowires (15, 22.5, and 30 nm). To remove the effect 
of the source and drain region, the phonon-limited mobility is 
calculated from the slope of the resistance [48]. The calculated 
phonon-limited mobility shown in Fig. 10(a) does not show a 
monotonic behavior with respect to diameter. It instead shows 
a peak at 2.5 nm and a decreasing trend as diameter increases 
to 3 nm. It can be explained by the argument that the higher 
subbands with a larger effective mass and scattering rates start 
to be populated with electrons as the diameter increases [52]. 

Fig. 10(b) shows the effective mobility calculated for the 
rough and perfect NWs compared to experimental results. 
To compare the simulation results with experimental results, 
we consider the fact that the ballistic mobility increases 
proportionally to the channel length and that the effective 
channel length is smaller than the gate length as described 
in Ref. [22]. The impact of IR on the effective mobility is 
still very small at diameter 4nm, but it increases as diameter 
decreases. 

The discrepancy between the experimental results [22] and 
the simulation results is noticeable even after phonon-limited 
mobility is included in the total effective mobility through 
Matthiessen's rule. We attribute this difference to the different 
electron density at which the effective mobility is calculated in 
the experiment and our simulation. In the experiment, supply 
voltage Vdd used is higher than that used in our simulation so 
that the effective mobility at the ON-state in the experiment 
is smaller than that in the simulation. In addition, the surface 
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quality of experimentally fabricated NWs may be different 
than that of the 2D plane surface which is adopted for our 
simulation. Nonetheless, IR scattering is crucial to understand 
the decreasing effective mobility in NWs of smaller diameter. 

Fig. 11(a) describes the importance of the quality of the 
interfaces on the electron mobility in NWFETs for different 
A m . The mean of the IR-limited mobility decreases with 
~ A~ 2 and it becomes comparable to the phonon limited 
mobility when A m = 0.2 nm. This indicates that poor surface 
quality can reduce the electron mobility significantly. 

When A m increases, not only does the effective mobility 
decrease on average as in Fig. 11(b), but also the standard 
deviation of the effective mobility increases. The standard 
deviations of the effective mobility with respect to the means 
of the effective mobility are given in percentages of 11.5 % 
(An=0.14 nm), 17.7 % (A m =0.2 nm), and 32.4 % (Z\ m =0.3 
nm). The standard deviations of the ON-current in Fig. 11(d) 
are, in percentages, 7.37 % (Z\ m =0.14 nm), 9.09 % (A m =0.2 
nm), and 21.6 % (zl m =0.3 nm) from the mean of ON-current. 
The variability of the ON-current due to IRS is smaller than 
that of the effective mobility. In conclusion, the quality of 
the interface between Si and Si02 is important in reducing 
the variability of the ON-current and the effective mobility in 
NWFETs. 

IV. Conclusion 

The effect of interface roughness between Si and Si02 on 
the performance of NWFETs has been studied through an 
atomistic 3D full-band simulation. The Si02 layer is modeled 
through VCA and its impact on the drain current is found to 
be significant. IRS is important in determining the threshold 
voltage of NWFETs. The experimentally observed [22] trend 
of nonlinear threshold voltage slope below 2.5nm is captured 
by the simulation. The significant reduction of the ON-current 
is observed in NWs with diameter 2 nm with more than 20 % 
reduction. 

Though the phonon scattering is found to be an important 
scattering mechanism in NWFETs, IRS cannot be ignored. 
IRS has a significant effect on both the ON-current and the 
effective mobility specially when the diameter decreases down 
to 2 nm. As the effect of IRS depends on the surface quality in 
NWFETs, the interface roughness scattering is a crucial factor 
to predict the device performance of NWFETs. 

Appendix A 

Calculation of Autocovariance function for Si0 2 

The random variables describing the rough Si02 surface, 
denoted as s(rQ, must be calculated from the Si atoms because 
they are inter-correlated with each other. Let us assume that 
one Si02 atom is connected to two Si atoms whose random 
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E p 


VsSCF 








SiQ 2 


-4.8 


1.83 


-2.27 


3.9265 


5.4655 


-0.3140 



TABLE I 

sp 3 Tight-binding parameters for the VCA S1O2 model. 



variables are s(ri) and s(r i+1 ). Then we approximate s(r^) 
by averaging the two random variables as 

s(r,;) + s(r i+1 ) 



(A) 



This approximation induces an error in the autocovariance 
of the Si0 2 atom (C^) as shown in the following calculation: 

= <*(ri)s(r$)) 



s(r-i) + s(n+i) 



s(rj) + s(r j+ i) 



^2 e -V2\ri-r } \/L m ^ + e V2Ar/L m + g -v / 2Zir/L m ^ 

1 + cosh (V2Ar/L m 



c. 



(5) 



where Eq. (1) and the relationship Ar = r,_|_i — r.; = rj + i —rj 
are used. Then, Eq. (5) can be further approximated with a 
Taylor expansion up to the second order as 

„2" 



C'ij 



a. 



1 + T 



(6) 



with x = y/2Ar/L m . If Ar w 0.1358 nm and L m = 0.7 nm, 
the difference between C,j and due to this approximation 

is x 2 /4 = Ar 2 /(2L 2 n ) « 1.88%. 

Appendix B 
S1O2 model and results 

The list of the sp 3 tight-binding parameters for Si02 is given 
in Table I and the bandstructure, including the imaginary band 
linking the lowest conduction band to the light-hole band at V 
in Fig. 12. The light-hole effective mass is strongly correlated 
to the electron effective mass in the sp 3 tight-binding model 
(TO;h=-0.3438mo). The heavy-hole mass does not affect the 
Si0 2 imaginary band linking the light-hole valence band and 
"lowest" conduction band at T. It is this latter imaginary band 
which largely determines the wavefunction decay constant in 
the oxide gap. Therefore, it is arbitrarily set to m/,/,=-1.5mo. 
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m =-1.5m 

hh 




-1 -0.5 0.5 1 

Im(k) (normalized) R e (k) (normalized) 



Fig. 12. Bandstructure of the fictitious VCA diamond-like Si02 material at 
k y =k z =Q. Real and complex bands are shown as well as the silicon conduction 
and valence band edges, the electron and heavy-hole effective masses of the 
oxide, and its band gap. 
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